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Abstract 

Chameleon fields are scalar fields whose mass depends on the ambient matter 
density. We investigate the effects of these fields on the growth of density perturba- 
tions on sub-galactic scales and the formation of the first dark matter halos. Density 
perturbations on comoving scales R < Ipc go non-linear and collapse to form struc- 
ture much earlier than in standard ACDM cosmology. The resulting mini-halos are 
hence more dense and resilient to disruption. We therefore expect (provided that the 
density perturbations on these scales have not been erased by damping processes) 
that the dark matter distribution on small scales would be more clumpy in chameleon 
cosmology than in the ACDM model. 
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1 Introduction 



Cosmological observations [1, 2] indicate that more than two thirds of the energy density 
of the Universe is in a component with negative pressure. Candidates for this missing 
energy, which is causing the Universe to accelerate, include a cosmological constant and 
scalar field models with equation of state w 7^ 1, often referred to as quintessence [3]. In 
these models, the scalar field is either decoupled from Cold Dark Matter (CDM) or couples 
to CDM but not to the baryons (coupled quintessence) 

In order to generate the present day acceleration the scalar field in these models must 
be evolving slowly and hence have a tiny mass, of order the present day Hubble constant, 
Hq ~ lO^^^eV. Since the mass of a quintessence field is very small, it can give rise to a new 
long-range force. Such a force has not been observed and consequently the coupling of the 
quintessence field to matter must be very small. Unfortunately, in effective theories derived 
from string theory nearly massless fields couple to matter with gravitational strength and 
would produce unacceptably large violations of the equivalence principle. Khoury and 
Weltman have proposed a scenario where a scalar field with gravitational strength coupling 
to matter can evolve on a Hubble timescale and generate the present day acceleration while 
evading all existing tests of gravity [4]. The key feature of this scenario is that the scalar 
field, which is dubbed the chameleon, has a mass which depends on the local background 
matter density. On Earth where the density is high, the Compton wavelength of the field 
is sufficiently small to satisfy all current tests of gravity, while on cosmological scales, 
where the density is tiny, the field has a much smaller mass and can drive the present 
day acceleration [5]. The field, however, is heavier than in standard quintessence models 
('^cham ^ H). In the solar system, where the density is many orders of magnitude smaller 
than on Earth, the chameleon is essentially a free field and mediates a long range force 
which could be detected by upcoming satellite experiments [6]. 

The cosmological history of the chameleon field was studied in Ref. [5]. It was shown 
that there is an attractor solution, analogous to the tracker solution in quintessence models, 
where the chameleon quickly settles into the minimum of its effective potential, and for a 
broad class of potentials and initial conditions the chameleon can satisfy all observational 
constraints. The evolution of the density perturbations was also investigated. It was 
shown that perturbations on scales smaller than the scale of the chameleon feel a larger 
effective Newton's constant which causes them to grow more rapidly. The length scale 
of the chameleon (0(100 pc) at present) is somewhat smaller than the scales probed by 
large scale structure observations. There has, however, recently been much interest in the 
properties of the first generation of dark matter structures to form in the Universe, and 
the possibility that they may leave an observable imprint in the present day dark matter 
distribution [7, 8, 9, 10, 11]. 

In this paper we examine the effects of the chameleon on the density perturbations 
on sub-galactic scales and the properties (formation epoch and over-density) of the first 
gravitationally bound structures to form. In section 2 we review the necessary aspects of 
the chameleon and its dynamics from Ref. [5]. In section 3 we extend the calculations of the 
evolution of density perturbations in Refs. [9, 11] to include the effects of the chameleon. 
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Finally in section 4 we examine the effect of the chameleon on the formation of small scale 
structure. 

2 The chameleon and its dynamics 

The action describing the chameleon field, %, matter and gravity has the general form 

where Mp\ = (SttG)"^/^ is the reduced Planck mass and Xm the various matter fields. 
The metrics governing the excitations of the matter fields are related to the Einstein 
frame metric g^i, via the conformal rescaling gjj^ = exp {2(3i(j)/Mpi)gni, where the (3i are 
dimensionless quantities of order unity [12]. Notice that the scalar field couples to all 
matter species including the baryons. As an example, this is the case for the radion field 
describing the interbrane distance in brane world models based on the Randall-Sundrum 
model [13] in which branes are nearby, for which l3i = jS = l/\/6. 

By varying the action with respect to x it can be shown [12, 14, 4, 5] that the dynamics 
of X are governed by the effective potential 

Kfr(x) = V{x) + exp (Ax/Mpi) , (2) 

where the matter density pi is defined as p\ = —g'^^Tji]^ exp (SPix/ Mpi) so that it is inde- 
pendent of X and conserved in the Einstein frame. If V{x) decreases monotonically with 
increasing x and /?i > 0, this potential has a minimum, Xmin, which increases with decreas- 
ing Pi. The mass of small fluctuations about Xmin increases with increasing Xmin so that 
the chameleon can evade local tests of the equivalence principle and flfth forces, due to the 
high local density. 

Fiducial potentials of the form 

y(x) = M4exp(MVx"), (3) 

with M — 10~^eV so as to produce the observed present day dark energy density (and also 
satisfy local tests of general relativity) were studied in Rcf. [5]. Assuming for simplicity a 
single matter component with density pm and coupling (3, the field value at the minimum 
of the effective potential satisfies 

/ M M p^exp(/3xmin(t)/Mpi) 

It was found that for a wide range of initial conditions the chameleon field reaches the at- 
tractor solution with x(^) = Xmin(^) before big bang nucleosynthesis and has a cosmological 
evolution in accordance with all observational constraints. 

The presence of a chameleon field can effect the growth of structure in particular on 
small scales. In the following, we will consider the modifications to the growth factor due 
to the gravitational effects induced by the chameleon. 
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3 Perturbation evolution 



We will follow Refs. [9, 11] and work in the longitudinal gauge, we will however use the 
notation of Ma and Bertschinger [15]. The perturbed line element reads 

ds^ = 0^(77) [- (1 + 2*) drf + (1 - 20) gijdx'dx'] . (5) 

The equations of motion for the CDM density contrast 5c and the divergence of the CDM 
velocity field ©c can be obtained from the energy-momentum conservation equation, which 
contains an additional term due to the exchange of energy with the chameleon 

7^(';r;. = %(5^x)T». (6) 

Here, i stands for the component i and T = T'^. The equations of motion are then given 
by (the dot represents the derivative with respect to r] and Ti = (da/dr;)/a) 

4 = -ee + 3(/) + /3(5x)-, (7) 
Oc = -(7i + /3x)©c + A;'(* + « . (8) 

The perturbed Klein-Gordon equation for the chameleon field % is given by 

/ d'^V\ ( dV \ 

[bxY + m{bx)- + W + a^-^ j bx + 2V1/ f — + /3pJ - mx = -^pAa^ , (9) 

and we will also need one of the components of the first-order perturbed Einstein equation 
(Poisson's equation in the sub-horizon limit) 

A;^ + 3n{<j) + n^) = ATtGa^dT^ . (10) 

From very early times onwards (before nucleosynthesis), the mass of the chameleon 
field is much greater than the Hubble expansion and the field sits in the minimum of the 
effective potential. Consequently the interaction scale of the chameleon field is always much 
smaller than the horizon H~^, and the evolution of perturbations on super-horizon scales 
is unaffected by the chameleon. Furthermore the chameleon docs not couple to radiation, 
since it is a traceless fluid, and the evolution of perturbations deep within the radiation 
dominated epoch is also as in standard cosmology [15, 11]. Once S^Pc 3> 5rPr (which for 
sub-galactic scales happens prior to matter-radiation equality), however, the dark matter 
terms dominate as the source in the Poisson equation, eq. (10), and the coupling of the 
chameleon to the matter density is now important. In particular, perturbations in matter 
will influence perturbations in the chameleon field and vice versa. 

On the sub-horizon scales we are interested in we can neglect the oscillations in the 
perturbations in the chameleon field, and x is also small as the field evolves along the 
minimum of the effective potential [5]. Following Ref. [11] we also neglect anisotropic 
stress, so that = ■0, and baryon anisotropics (but not the baryon density). The latter 
assumption is valid on small scales (/c > /cb ~ 10~^pc~^) for 2; > 2^5 ~ 150 as prior 
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to this residual electrons allow transfer of energy between the photon and baryon fluids 

and thermal pressure prevents the baryon perturbations from growing [16]. Following 
Ref. [5, 11], eqs. (7-10) can be combined to give the following equation for the evolution of 
the cold dark matter density contrast: 



5c = -m, + - 



3 p, 



■c 



2 pc + Pj 



1+ 



1 + ^ 



4, (11) 



with V" = d'^V/d'x^. The effects of the chameleon manifest themselves in the second term 
in the square brackets. The chameleon field operates on length scales smaller than 

Acham(i) = -^Jyjj , (12) 

or equivalently for comoving wavenumbers larger than 

fccham(t) = " ■ (13) 
^chaml,''^ 

For k <S /ccham(^) the terms in the square bracket in eq. (11) above are well approximated 

by 1, and the CDM density contrast evolves as in standard cosmology. For k ^ kchaxait), 
they are well approximated by (1 + 2f3^) i.e. on these scales the growth of perturbations is 
governed by an effective gravitational constant given by G(l + 2/3^). 

For k 3> /Ccham, Eq- (11) can be re- written in terms of |/ = a/ceq to give 

y(y + 1)5': + (l + Iv) = ^ (l + 2/9^) (1 - h) <^c, (14) 

where ' = d/dy and we have introduced the baryon fraction = fib/f^m- The solution to 
this equation is a superposition of Legendre functions of first and second kind of order u: 



5c{k, y) = B,{k)P, {Jl + y)+ B2{k)Q, \/l + 1/ • (15) 



This is the same as in the ACDM case [9, 11], but the degree u of the Legendre functions 
in the standard case is given by 



-l + ^l + 24(l-/b) 

J^GR = ^ , (16) 

whereas we have, for k ^ A;cham(^), 

-1 + ^1 + 24(1 + 2/32) (1-/,) 

i^cham — ^ ■ [i-l) 

For the best fit WMAP ACDM model fb = 0.17 [1] so that z/qr = 1.8, whereas for /? = 1 

i^cham = 3.4. For z < Zi, 150 the baryons follow the CDM (which is equivalent to setting 
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fb — in eqs. (14)-(17)) and i/qr = 2 for the standard cosmology and i^cham = 3.8 for the 
chameleon with (3 — 1 {as found in Ref. [5]). For y ^ 1, equivalently z <^ Zeq, Sdy) grows 
as a*^/^. The increase in the growth rate due to the chameleon is large compared with the 
suppression in growth due to baryons which occurs for < z < z^q. We therefore now 
neglect the baryons and set /b = in eqns. (16) and (17) above. 
The full asymptotic late time solution is [9, 11] 



(18) 



where Co is the superhorizon limit of the curvature perturbation on uniform density hyper- 
surfaces and the constants c(i/) and b{i') are found by matching the early time radiation 
domination {y <^ 1) expansion of eq. (15) to the sub-horizon limit of the general radiation 
domination solution [9, 11] 

r(i + 2z/) , , 

and 

where T'^u) is the derivative of r(z/) with respect to u. For i/QR — 2, c — 1.5 and b — —1.7, 
while for Ucham = 3.8, c = 3.9 and b = —2.8. 

We now calculate the evolution of kchamit) with time, and consequently the scales on 
which the chameleon effects the growth of the density contrast. During matter domination 
V{xm.in) ~ const and Pmexp (/9xmin/Afpi) oc a'^ so that, using eq. (4), 

Xmin(t)=Xmin(to)(l + ;^)-'/^"^'^ • (21) 

At present PmCxp (/30min/-^pi) ~ ^(^min) ~ so that [5] 

/ X !/(-+!) . ^ .n/(n+l) 
Xmin(io)=U (t:^ Mpi, (22) 



\PJ VMpi/ 
Using eqs. (12) and (13) the scale of the chameleon field varies as 

Acham(i) = Aeham(io)(l + ^)^-'^"+'^/'("+'^ (23) 

,'1 I \(n+4)/2{n+l) 

^cham(i) = A;eham(to)(l + ^)^"+^^/^^"+^^- ^ J , ■ (24) 



where 



. 1 ^MpiV"+^)/^(-+^) 



The characteristic wavenumber A;cham(^) decreases with decreasing red-shift. At z, 



eq 



(1 + ^eq)^''"^''^/^^""^^^ 
/f^cham(^eq) = 7 Tr\ ' (26) 

'^chaml.roj 
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where (l + ^eq) = 24000Jli„/2'^ ~ 3700. For fiducial parameters n = /3 = 1 and M = lO'^eV, 

Acham(to) = 250 pc and Acchaml^cq) = 120 pc"^ 

For scales k < A;cham(-2cq) the density contrast growth law is only modified once k > 
kcha.m{t). We denote the redshift at which this happens by -Zmod, which is given by 

(1 + ^mod) ^ (Acham(io)^)'^"+'^^^"+'^ ■ (27) 

For z <C -^eq density perturbations grow as 5 oc a''^^, where for standard general relativity 

and k ^ fccham(^)i ^ = ^gr while for k ^ fccham(^)i ^ = ^cham- The, scale dependent, 
additional growth in the linear density perturbation due to the chameleon is given by 

where z^od ~ -^eq for k > A;cham(^eq) and is given by eq. (27) for k < /ccham(^eq)- 

In CDM cosmologies structure forms hierarchically (large halos form via the merger and 
accretion of smaller subhalos) and at least some substructure is expected to survive. There 
must be some cut-off in this process however; if the density perturbation spectrum extended 
down to infinitely small scales the contribution of density perturbations to the local energy 
density would diverge [7]. For weakly interacting massive particles (WIMPs) damping 
processes [17, 7, 8, 9, 11, 18], namely coUisional damping (due to clastic interactions with 
radiation) and free-streaming, produce a cut-off in the (processed) power spectrum at 
^cut ~ lpc~^ (i.e. fluctuations on scales k > kcut are erased). Ref. [11] found a range of 
values kcut ~ 0.4 — 4 pc~^ for benchmark models spanning the range of plausible WIMP 
properties. Dirac hke WIMPs, where elastic scattering is mediated by exchange, have 
values of fccut at the lower end of the range. Majorana WIMPs for which Z° exchange is 
suppressed, for instance neutralinos, have a wider range of kcut values with more massive 
WIMPs having larger k^ut- 

We plot A;cham(^) and Zraod{k) for n = 1 and n — 2 and also the range of /Ccut values in 
figs. 1 and 2 respectively. For n — 1 kcham.{zeq) is larger than the upper end of the range of 
kcut values and the evolution of the surviving perturbations is initially unaffected. kcha.m{z) 
decreases sufficiently rapidly with decreasing red-shift, however, that -2mod(^"cut) ^ for 
the entire range of kcut values and the growth law of small (physical) scales is modified. 
For n — 2 /ccham(-2eq) IS SO large that the chameleon scale is beyond the cut-off scale even 
at late times and the growth of surviving perturbations is completely unaffected by the 
chameleon. We should emphasise, however, that WIMPs arc not the only viable CDM 
candidate. There are a large number of other candidates [19] (including the arguably 
equally well motivated axion) whose microphysics have not yet been studied and which 
may have substantially different cut-off scales. 

4 Small scale structure 

The enhancement of the growth rate of the CDM density contrast on small scales due 
to the chameleon means that these scales will go non-linear, 5c ~ and collapse to 
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Figure 1: The red-shift dependence of the characteristic wavenumber of the chameleon, 
^cham(-2), for /5 = 1, M = 10~^eV and n = 1 (sohd hne) and n = 2 (dotted Une). The 
dashed hnes show the range of values of kcut for plausible WIMP models from Ref. [11]. 

form structure earlier than in standard cosmology. The red-shift at which the first typical, 
l-cr, fluctuations on comoving scale R go non-linear -2ni(-R) is defined via (T(i?, ^ni) = 1, 
where a{R, z) is the mass variance calculated by integrating the density perturbation power 
spectrum multiplied by the Fourier transform of a window function. This calculation is not 
possible in this case as we only have the evolution of density perturbations on sub-galactic 
scales, and the change in the growth law due to the chameleon means that we can not use 
the value of erg measured by WMAP to evade this problem, as in Ref. [9, 11]. We instead 
estimate the effect of the chameleon on the red-shift at which scale go non-linear by using 
the approximations (5cham(^, -2ni,cham(^)) = 1 and (5gr(/c, -Zni,GR(^)) — 1 and using the values 
of -2ni,GR(-R) calculated in Ref. [11] with the approximation k ~ 1/R. Using eq. (28) we 
find 

l + ^nl,cham(A;) ^ / l + ^inod(fc) \ ^'-^"^ '"^"^^^^ 
l + ^nl,GR(A;) Vl + ^nl,GR(A;)y 

The resulting values of 2;ni^cham(^) are plotted in Fig. 4 For scales with k < A;cham(^eq) we 
assume that v, and hence the growth law for 5c, changes abruptly at ^^mod- The change 
would in fact occur smoothly as the chameleon term in square brackets in eq. (11) increases 
smoothly from 1 to (1 + 2/3^). This approximation is reasonable however given the other 
uncertainties involved in the calculation. 
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Figure 2: The red-shift at which the growth law of density is modified due to the chameleon, 
^mod(^); for /3 = 1, M = 10~^eV and n — 1 (solid hne) and n = 2 (dotted line). The short- 
dashed lines show the plausible range of kcut values for WIMPs. 

A scale dependent primordial power spectrum, with spectral index > 1, as produced 
by, for instance, false vacuum dominated hybrid inflation would also result in a larger than 
standard density contrast on small physical scales and hence earlier structure formation. 
We therefore also plot -Zni.oR as calculated in Ref. [11] for a false vacuum dominated hybrid 
inflation model which produces a primordial power spectrum with Ug = 1.036, which is 
the maximum scale dependence allowed by the WMAP and 2dF data [20]^. For n = 2 
^ni,cham(^) = z-ai^GR{k) foT all k < lO^pc"^. For n = 1, foi k > 0(lpc~^) the growth law 
is modified sufficiently early that 2;ni,cham(/i;) ^ -2ni,GR(^)- The rapid increase of 2;ni,cham(/i;) 
with increasing k is caused by the large change in the index of the growth law by the 
chameleon, 5c oc a^'^ compared with (5c oc a for standard cosmology. This resulting change 
in 5c on small physical scales at late times is far larger than that produced by a primordial 
power spectrum with rig = 1.036. 

The physical properties of the first generation of WIMP halos were estimated in Ref. [11] 
using the spherical collapse model. The physical size of the halos post collapse is propor- 
tional to (1 -I- Zn\)~^ and hence the density contrast is proportional to (1 -|- Zni)^ i.e. the 

possible interaction between the chameleon and the inflaton field docs not change the result for Us, 
since the chameleon field is rather heavy during infiation and has no influence on the effective mass of the 
inflaton field [21]. 
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Figure 3: The red-shift at which typical 1-a fluctuations collapse to form structure, 
-2^ni,cham(^), foi /5 = 1, M = 10~^eV and n = 1 (sohd line). The dotted line is 
^ni,GR(-R = 1/^) and the short-dashed lines show the plausible range of fccut values for 
WIMPs. The long-dashed line shows Zni,GR{R = for a primordial density perturba- 
tion power spectrum with spectral index Ug = 1.036. 

earlier halos form the more over-dense they are at later times (reflecting the higher matter 
density at the time they form). In fig. 4 we plot the present day overdensity corresponding 
to typical halos, which form from 1-a fiuctuations. A, 

_ 2MiR) 

- f r(i?)3p^(to) ' ^'"^ 

as a function of comoving scale where M{R) = 1.6x lO^M0(f2jn/i^/O.14)(i?/pc)3 is the mean 
mass within a sphere of comoving radius R and r{R) = 0.53i?/(l -|- ^ni) is the physical, 
post collapse, radius of a halo which forms from a typical fiuctuation with comoving size 
R at red-shift ^ni- For the chameleon we make the approximation i? ~ l//c. The first halos 
to form in chameleon cosmology are significantly more concentrated than in standard 
cosmology (provided A'cut > 0{lpc~^) if the CDM is in the form of WIMPs) and are 
hence more likely to resist disruption by dynamical processes (such as tidal disruption and 
interactions with stars) . We therefore expect that the present day dark matter distribution 
on small scales would be more clumped than in standard cosmology and this could be 
detectable via axion detectors or WIMP direct and indirect detection experiments (the 
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survival probability of the first dark matter halos even in standard cosmology is the subject 
of ongoing studies [8, 10, 22]). This alone would not provide a 'smoking gun' for the 
chameleon, as other modifications of general relativity could also lead to enhanced growth 
of small scale density perturbations and hence small scale structure. However, combined 
with a detection of modified gravity within the solar system [6] the present day density 
distribution could be used to probe modifications of gravity, such as those due to the 
chameleon. 




-2-10 1 2 

logjo (R /Ipc) 

Figure 4: The present day overdensity corresponding to typical fluctuations. A, for the 
chameleon with P — 1, M — 10~^eV and n — 1 (solid hue) and for standard cosmology 
with primordial power spectra which are scale independent (dotted) and Us = 1.036 (long- 
dashed). The short-dashed hues show the range of l//ccut values for plausible WIMPs. 
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